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Abstract 

s-channel discontinuity of quark-quark scattering amplitude with gluon 
quantum numbers in the t channel and negative signature is calculated in 
the Regge kinematical region in the two-loop approximation. Using this dis- 
continuity and assuming that the Regge asymptotic behaviour is given by 
the Reggeized gluon contribution, we calculate the gluon trajectory in the 
two-loop approximation. Remarkable cancellations lead to the independence 
of the trajectory on properties of the scattered quarks, confirming the gluon 
Reggeization. 
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1. INTRODUCTION 

Perturbative QCD is widely used nowadays for the description of semihard | I 
as well as hard processes ||. However, whereas for the later ones the theory is well 
developed and understood, for the former ones we have a lot of unsolved theoretical 
problems. The applicability of perturbation theory improved by the renormalization 
group for a hard process having a large typical virtuality Q 2 is justified by the 
smallness of the strong coupling constant a s (Q 2 ). Contrary, for semihard processes 
another essential parameter appears: the ratio of the typical virtuality Q 2 to the 
square of the c.m.s energy s of colliding particles. At sufficiently high energy the 
value x = Q 2 /s becomes so small that it is necessary to sum up terms of the type 
a"(/n(l/x)) m , with m < n (for the scattering channel which is considered here). Up 
to now this problem is solved || in the leading logarithmic approximation (LLA) 
only, which means summation of the terms with m = n. 

The results of LLA have two serious disadvantages. Firstly, the Froissart bound 
atot < c In 2 s is violated in LLA. In fact, calculated in LLA, the total cross section 
a tot A grows at large c.m.s. energies as a power of s: 

-LLA 6 ri\ 
<Jj , ~ 



tot r, j 

V in s 



where, for the gauge group SU(N) (N = 3 for QCD), with gauge coupling constant 
9 («. = 

oj = i-N\n2. (2) 

In terms of structure functions this means their strong power increase in the small 
x region. The Froissart bound is violated in LLA because the s-channel unitarity 
constraints for scattering amplitudes are not completely fulfilled in this approxi- 
mation. The problem of unitarization of LLA results is extremely important from 
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a theoretical point of view. It is concerned in a lot of papers (see, for example, 
Ref. 1). 

Another disadvantage seems to be even more important from a practical point 
of view, since the results of LLA are applied to the small x phenomenology (see, 
for instance, Ref. M). There is an uncertainty of the argument of the running 
coupling constant which appears because the scale dependence of a s is beyond of 
the accuracy of LLA. This uncertainty diminishes the predictive power of LLA, 
permitting to change strongly numerical results by changing a scale. 

Therefore, the problem of the calculation of radiative corrections to LLA be- 
comes very important now, as it gives us the possibility to fix the scale dependence 
of the coupling constant, to reduce the uncertainty of the predictions of LLA and 
to determine a region of its applicability. 

A solution of this problem can be strongly simplified || by using the Reggeiza- 
tion property of the non-Abelian SU(N) gauge theories. It was proved || [!]] in LLA 
that gauge bosons are Reggeized in these theories with trajectory 

j{t) = l + tj{t) , (3) 

where in the leading approximation 

m/ \ g 2 t N r d D ~ 2 k ± 
u(t) = u^it) = - / — j^r • (4) 



(2tt) (D - 1) 2 J kl{q-k)[ 
Here q is the momentum transfer, t = q 2 »s q±, and D = 4 + e is the space- 
time dimension. A non zero e is introduced to regularize Feynman integrals. The 
integration in Eq.(^) is performed over the (D— 2)-dimensional momenta orthogonal 
to the initial particle momentum plane. 
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The results of LLA are summarized || in the Bethe-Salpeter type equation for 
the t-channel partial amplitudes. The problem of the calculation of corrections to 
LLA, therefore, can be set up as the problem of the calculation of corrections to the 
kernel of this equation ||. The kernel is expressed in terms of the gluon trajectory 
and the Reggeon-Reggeon-gluon (RRG) vertex. The corrections to the vertex are 
known || |9|, therefore the calculation of the contribution u^(t) to the trajectory 
in the next (two-loop) approximation has become most urgent. 

In this paper we present results and details of the calculation of the two-loop 
correction to the trajectory for the real case of QCD with massive quark flavours. 
The result for uj^ 2 \t) in the massless quark case was obtained earlier |T(|. The paper 



is organized as follows. In Sec. II we discuss the method of calculation. In Sec. Ill 
we calculate the contribution of the two-particle intermediate state to the s-channel 
discontinuity of the quark-quark scattering amplitude. An analogous calculation is 
performed in Sec. IV for the contribution of the three particle intermediate state. 
The final expression for the correction uj^it) is obtained and discussed in Sec. V. 

2. METHOD OF CALCULATION 

The method is based on using s-channel unitarity. The two-loop contribution 
to the gluon trajectory can be obtained from the s-channel discontinuity of an 
elastic scattering amplitude with gluon quantum numbers and negative signature 
in the t channel calculated in the two-loop approximation with the accuracy up 
to a constant. Indeed, let us consider such an amplitude for a process of the type 
A + B — ► A' + B' at large s and fixed t. Assuming the gluon Reggeization, the 
amplitude takes the factorized form 



, \\ A' B' . S 

"^8 ) — 1 A' A 

/ AB AA t 



-t) + b 



^B'B 5 (5) 



where T l A , A are the particle-particle-Reggeon (PPR) vertices. They can be written 

as 

Ta'A = 9(A'\T\A)(T^ a + ^a) , (6) 

where (A'|T l |A) stands for a matrix element of the colour group generator in the 
corresponding representation (i.e. fundamental for quarks and adjoint for gluons), 
T A ) A and T A ) A are respectively the Born and the one-loop contributions to the 
vertices. Using Eq.(H) and the decomposition 



one can present the two-loop contribution to the s-channel discontinuity 
in the form 



(7) 



A B 
AB 



A 8 { >) ab (two-loop) 



2nis . 



)g 2 (A'\T\A){B'\T\B) 



r%(^(t)r in (4)r^ + (r&rgk + r^) B )^\t) + rg^Wflrg? 



-t 



B 



Since the one-loop corrections T A ) A to the PPR vertices became available P, |TT], [12|, 
the only unknown quantity in the RHS of Eq.@ is the two-loop contribution uj( 2 >(t) 
to the gluon trajectory. Consequently, one may obtain it from the expression of the 

A'B' 



discontinuity 



(^8 h) ) AB (tWO-loop) 



calculated in the two-loop approximation 



with accuracy up to a constant. The calculation of this discontinuity is the main 
content of this paper. 

By definition the trajectory should not depend on a particular type of scattered 
particles, therefore we have a freedom in choosing these particles. In this paper we 
consider the process of the quark-quark scattering for calculating the correction to 
the trajectory. 



In LLA helicities of each of the scattered particles are conserved, so that in the 
helicity basis we have 

r®4 = Sx Al x A • (9) 

On the contrary, in higher orders PPR vertex T^, A may contain another spin struc- 
ture. Because of the parity conservation, the one-loop correction T^) A can be written 



as p| 



rSi* = 5a Ai a/2W + Sx A ,-x A ,r ( A l(t) (io) 

if the relative phases of the states with opposite helicity are appropriately chosen. 
Both quantities T^{t) are calculated in Ref. |12| . For our purpose the knowledge 
of r^(f) is sufficient. In fact, the correction u^(t) in Eq.(H) is multiplied by the 
helicity conserving vertices, then for its determination only the helicity conserving 
parts of the amplitudes entering Eq.(|) are necessary. 

Let us write the helicity conserving part of the discontinuity of the amplitude 
in the LHS of Eq. (||) in the following form: 



, s\A'B' 

As Ub ( tw °- lo °P) 



™ = 9 2 (A'\T\A)(B'\T\B) (-^) A s , (11) 

where the superscript (+) in the LHS means helicity conserving part. Then Eq.(H), 
together with Eqs.(|9|) and QlCf), gives us 

J*\t) = A 8 — (c«(t)) 2 in (^) - (r«(t) + r«(t)) . (12) 

The only unknown term in the RHS of Eq.(^) is the discontinuity A s . We calculate 
it below by using the s-channel unitarity condition. In the two-loop approximation 
only two- and three-particle intermediate states do contribute. This allows to divide 
the discontinuity under consideration into two parts respectively: 

A s = Af + A® . (13) 



3. TWO-PARTICLE CONTRIBUTION TO THE DISCONTINUITY 

Let us consider the two-particle contribution in the unitarity relation: 



aA'B' 
^AB 



= I 

5(2) 



/ d$ 2 { PA + PB ; PAl ,p Bl ) £ A% Bl A A % Bl . (14) 

Since we are considering the quark-quark scattering amplitude, the intermediate 
particles can be only quarks and the summation is performed over their spin and 
colour states; <i$ 2 is the two-body phase space element. Its general expression reads 



(n \ n 
p-Y,pi n 
i=i ) i=i 



(15) 



In order to obtain the two-loop contribution to the discontinuity we need to take 
one of the two amplitudes in the RHS of Eq. (|14]) in the Born approximation and 
the other one in the one- loop approximation (see Fig. 1). Both amplitudes should 
be calculated in the region of asymptotically large c.m.s. energies y/s and fixed 
momentum transfers. In this region the amplitudes acquire a very simple structure 
in the Born approximation: 



^(Born) = (A 1 \r\A)(B 1 \T]B)5 XA , XAi 5 XB , XBi ^ , (16) 

where t\ = (p A — PaJ 2 - The amplitude A^i^ may be obtained from Eq.([T6|) by 
evident substitutions. 

The one-loop contributions to the amplitudes A AB Bl and A^^i have a much 
more complicated form. Fortunately, it is sufficient to take into account only those 
parts of these contributions which conserve the helicity of each particle and have 
negative signature in t\ and t[ channels correspondingly (t[ = (p A > — PaJ 2 )- Indeed, 
let us consider the spin structure. We are interested in the helicity conserving part 
of the discontinuity, so that only such a part should be kept in the product of the 
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2g 2 s 



Born and one- loop contributions in Eq.(14). Since in the Born approximation the 



helicities are conserved, as Eq. flTBI) shows, the only helicity conserving part must 
be taken in the one-loop contribution. Now let us turn to the colour structure and 
signature. In the Born approximation (16) the amplitudes evidently have gluon 



quantum numbers in t\ and t<i channels and negative signature (in the region of 
asymptotically large s it simply means that they are odd functions of s). The one- 
loop contribution to the amplitudes A AB Bl an d A^l^ , besides such a part, has also 
a positive signature part, with colour singlet as well as octet state in the t\ and t[ 



channels. This part, however, is purely imaginary [12|, whereas the Born amplitude 



is real; therefore, their product cancel in the unitarity equation (|T~[|). 

(2) 

As for the result, to calculate the discontinuity A^ we may use, for the ampli- 
tudes A AB Bl and A A }& in Eq.(|14|), the representation (|5|) with the trajectory oj(t) 
and vertices T C A , A given by Eq.((4|) and Eqs. (|]), (^) and fllCf ) respectively. In the 
one-loop approximation we get 

^(one-loop) = {A 1 \T\A){B l \T i \B)5 XA ^ Ai 5 XB , XBi 

x & {u^fa) [in ( JL) + In (ZLJ\ + 2 [i^fe) + r«(ti)] } + •••, (17) 

where dots denote terms which do not contribute to the discontinuity A^ 2 ). The 
amplitude A A i^ (one- loop) is obtained from Eq. (p~7f ) by evident substitutions. Using 
the Born (|T6| ) and one-loop (|1^) amplitudes and summing up in Eq.([L4|) over spin 
and colour states of the particles A\ and B\ leads to the spin and colour factor 

(A'\T^'\A)(B'\T i T'\B)8 XA * A ,8 XBtXB ,. 

Let us remind that we need to calculate not the discontinuity (|14D itself, but the 
corresponding discontinuity for the amplitude with the octet colour state in the 



t channel and negative signature. The calculation of the discontinuity for this 
amplitude requires antisymmetrization with respect to the change A — ► A', A' — ► A 
or, equivalently, B — > 5', 5' — > 5, that in our case reduces to the change of the 
order of the group generators into the factor (A'\T l Ti\A) . Using the relations 

N . 



rpl rpj 



if 



k 1 



f.. k T l T j = -—T k , 



;is) 



where the coefficients / ij are the group structure constants, the colour factor be- 
comes 

N 

- — (A'\T\A)(B'\T\B) . 

Notice that the negative signature automatically leads here to the colour octet state 
in the t channel. Using the two-body phase space element 

l S D ~ 2 ki± 



with qi = pa± — Pa, from Eqs.flHD, (|16D and (|17|) we get 



(19) 



{ As Jab ( tw °- lo °P) 

d^qii. 



(+) 



(Qi ~ <?)i<7i 2 ± 



5(2) 

uU(q! ± )ln 



g\A'\T l \A){B'\T l \B) 



2nis\ g 2 Nt 



t ) (2tt) d - 1 



, ,+ri + i(^) + rW(^; 



(20) 



Comparing Eq.(|20|) with Eqs.(|TlD and (|T3|), we obtain 



9 2 Nt 



Qi± 



XL) 



(21) 



(2vr)^ J ql x ( gi - q)l 
The helicity conserving part of the one-loop correction to the quark- quark- Reggeon 

vertex r^lit) is calculated in Ref. [l!|. For our purpose it is convenient to express 

it as a sum of three parts having different flavour-colour dependence: 



TqqW = a f(t) + oq(*, m o) + a g (t, ™q) 



(22) 



Here the first term 



a f (t) 



(4tt) 



D 



2,-,o 



x(l — x) 



m 2 — x(l — x)t 



2-4f 



(23) 



is the contribution of the quark loop, the summation being over the quark flavours. 
The second term 



a Q {t,m\ 



g 2 1 / D 
z r 2 

(47r)T2AT V 2 




dx 



D-A 



2m 2 



+ 



2 \4r-2 



(24) 



,D-3 4 7 (-D-3) 
is connected with the vertex and the quark self-energy diagrams. Eq. (^4j) is obtained 
from Eqs.(39) and (68) of Ref. [12| with the help of the identity 



(D-3) / dx 



tx(l — x) 



[m 2 — x(l — x)t] 



3-T 4 7o 



1 /-i , t(D-4)+4m 2 
' ctx- 



[m 2 -x)tf 2 



(m 2 )' 



which in turn follows from the evident identity 



dx— 



x 



(m 2 )' 



(25) 



(26) 



/o cte \j m 2 _ x (i _ x ) t y 
Finally, the last term, proportional to N, comes from the two gluon exchange and 
quark self-energy diagrams; it can be written as 



Here a g (t, 0) is given by Eq.(63) of Ref. @, 

2jV r (2-f)r 2 (f -l) r 

a, f , 0) = D 2 ^ V Z P - 3 

(4tt)- (-t) ~ T (D - 2) I 



^3 



D 



-2^(--2)+^(l) 
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4(D-1) D-4 4 



(27) 



f28) 



where ip(x) is the logarithmic derivative of the gamma function: 

it \ T ' ix) 
ip(x) 



T(x) 



In turn, the second term 5 S (£, m^) comes from Eq.(64) and from the part propor- 
tional to N in Eq.(67) of Ref. |12| . It is convenient to modify its form using the 
identity 



J J dx\dx 2 Q{\ — X\ — x 2 ) ^ 



8 -2 8 



3x 2 dx X ) [ m 2 x 2 _ (1 _ Xl _ X 2 )f-^ 



2\¥—2 



D-2 



m 2 



which provides us the relation 



2 / / dx\dx 2 
Jo Jo 



i r i 



tx 2 (l - xi - x 2 )8(l - xi - x 2 ) 2 2 «,2--i 

tt H m p 

[m 2 x{ — tx 2 (1 — Xi — x 2 )] 2 — 

0(1 - xi - x 2 ) [(§ - 2) tei(l - x x ) + 2{D - 3)m 2 xj 



dxidx 2 ^ 

/0 Jo \m 2 x\ — tx 2 (1 — Xi — x 2 )] 2 

Using this relation one arrives at 



5g(t,m 2 ) = 9 N n T !' 3 - 



(47T) 



y) I ^ / dxi( ^ x2 ^^ ~ Xi ~ X2 ^ 



t(l -x x ) 



1 — Xi + 



D-2 



( 



-x 



3— — 

[-tx 2 (1 - Xi - x 2 )j 2 J 



m Q x i ~ t x 2 (1 ~ x l ~ X 2) 

2m 2 Q xi 



1TLqX 2 — tx 2 (1 — Xi — x 2 ) 



6 2 



2{m 



2^-2 



(D-3)(D-A) 



(29) 



(30) 



(31) 
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4. THREE-PARTICLE CONTRIBUTION TO THE DISCONTINU- 
ITY 

Let us now turn to the calculation of the contribution coming from the three- 
particle intermediate state in the s-channel unitarity condition. Since we are consid- 
ering the discontinuity of the quark-quark scattering amplitude, in the intermediate 
state we face the same quarks with an additional gluon. The unitarity condition in 
this case reads 

Aii'] sm =iJd^(p A +p B ;p Al ,p G ,p Bl ) £ A% G ^Af B ? Bl ■ (32) 

where G is the produced gluon. For calculating the discontinuity it is convenient 
to use Sudakov variables. Let us introduce light-like momenta p- and p + close to 
the initial momenta pa and ps respectively: 

m 2 A m 2 R 

p A =P-^ P+ , p B =P+^ P- , 

s s 

p 2 _ = p 2 + = , %P+P- = s . (33) 
Thus any vector pi can be decomposed as 

Pi = fop- + oiip + + pix , (34) 

so that 

p 2 = sai(3i+p 2 i± . 
In the region under consideration, i.e. for 

s 3> \t\ ~ m 2 A ~ m 2 B , (35) 

the essential kinematics in Eq.(^) is defined by three conditions: 

*) \Pm±\ ~ IPbuI ~ \Pg±\ ~ 1*1 > ( 36 ) 
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which means that all transverse momenta are limited; 

ii) f3 Al ~ a Bl ~ 1 , ai Al ~ (3 Bl ~ -j , (37) 

which states that the order of magnitude of quark longitudinal momenta does not 
change; 

in) - < Pa < 1 , - < olg £ 1 , (38) 

s s 

according to which the range of variation of the gluon rapidity is large. Since 

sa G (3 G = -p 2 G± ~ \t\ , 

only one independent variable varies strongly in the multi-Regge region, leading to 
a contribution proportional to ln(s) 0. 

In LLA, where the discontinuity should be calculated with logarithmic accuracy, 
only the multi-Regge kinematics could contribute so that the region was defined by 
Eq.© and by 

S 

— < fa < 1 , — < oi G < 1 ■ (39) 

s s 

Now we need to keep a constant term together with the logarithmic one, therefore, 

besides the contributions of the region ( p8|) we must also calculate those ones of the 

fragmentation region of the particles A, A' and B, B' . Let us start with the first of 

them. The region is determined, together with Eq.(^), by the relations 

Pax ~ 1 , a Bl Ril, a Ax ~ Pbx ~ — 5 

s 

/3 G ~ 1 , a G ~ — . (40) 
s 
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Within the accuracy here required, the gauge invariant expression for the amplitude 
•^ab Bi m this region is given by 



Ajg* = ^u(p Al ) {(A^T^A)^ ^ -^ G + mA) 



-2(p A p G ) 



G 



+ (Ai\T c T j \A) f G + ^ G + mA) jf + - ifij C (Ai \T\A) 
2{Pa 1 Pg) 

x {pAi 2 _ pA)2 [fc(P + Pc) - MP + eh) ~ A ( [pa - PM + Pb) ^ 

x u( PA )e* G c (B 1 \T^\B)- ^u(jp Bi y_u{PB) ■ (41) 

All Feynman diagrams for the gluon emission in the quark-quark scattering (see Fig. 
2) contribute to the amplitude ( ffip and one may easily recognize the contribution 
of each diagram in the Feynman gauge. Let us remind that, if a diagram contains 
a gluon line connecting two parts with strongly different Sudakov variables, then 
for calculating its asymptotic contribution it is convenient to decompose g^ u in the 
gluon propagator in the form 

f = 2 (MlM + 9f . (42 ) 



The asymptotic contribution is given only by the first term in the RHS of Eq. ( fl2| ) . 
supposing that the Sudakov variables of /3-type (a-type) at the /i-vertex are much 
smaller (much larger) than those at the v-vertex. Clearly, this trick is used for the 
gluon lines with momentum q 2 in the diagrams of Fig. 2(a,b,c). It is also used 
for the gluon lines with momentum qi in the diagrams of Fig. 2(d,e), because the 
dependence on the gluon momentum is factorized for these diagrams, in the same 
way as for the soft gluon emission, and after this we have all the conditions for 
applying the trick. 
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We could perform all calculations in a gauge invariant manner, but it occurs 
that a suitable choice of the gauge makes our task easier. This happens in the axial 
gauge 

e G P+ = • (43) 

Using then this gauge condition, let us express the amplitude (ffl|) in terms of the 
Sudakov variables. First of all we introduce the simpler notations 

k = p G , e = e G , q 1 = p A - p Al , q 2 = p Bl ~ Pb , 

k — qi - q2 — ap+ + @P- + k± . (44) 

Moreover, by means of the gauge condition ( f4"3"D and of the transversality require- 
ment ek = 0, we can express the polarization vector in terms of its transverse 
part: 

e = e±-p + -. (45) 

P+k 

Hence all the invariants necessary to calculate the amplitude A^^ Bx acquire the 
following expression in terms of the Sudakov variables: 

2 P+ k = s(3 , 2p„k = -S , 2p +PAl = s (l - (3) , 
2p. PAl = ^4±, 2p A k- m ^ 2 - k - 



1-/3 ' " (3 
_ - (k + (3q 2 )l _ ~ m\(3 2 , _ , 

*PM*- (3{1 _ (3) , ii- 1 _ f3 . <? 2 -g 2 x- (4b) 

Using now the polarization vector (£E|), the invariants fl4"6|) and the commutation 
relation fll8|) for the term of Eq . (f4l"|) containing (y4 1 |T c T- J |A), through a simple 
calculation of spinor algebra, we arrive at 

A A A f B ' = 2g*u{p Al ) {(A 1 |T^T C |A)[L(^) - L(k ± + f3q 2S _)} 
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HMA|T*|A>[L( gi± )-L(A ;± + ^ . (47) 

Here we have put 

L(k ± ) = L,(k ± )eT , L»(k ± ) = 7 " (m ^~/_^ + 2 ^ • (48) 



In order to go from Eq.([y]) to Eq.fl47D we must take into account that in the helicity 
basis one has 

u(pbi)i/-u(pb) = s6 Xb ,x Bi . (49) 

The corresponding expression for the amplitude A A l^, Bl can be obtained in the 
same way and reads 

A A A }^ = 2g 3 u(p Al ) {(A x \T^T c \A'){L{k L + (3q ± ) - L(k ± + (3q 2A _)) 

-if ijc (A 1 \T i \A')(L(q 1± - (1 - (3)q ± ) - L(k ± + (3q 2i _))} 

x tf + u{p A ,)e*S 1 m^Wi^ . (50) 

[Q2 — g)± 

Let us note that this expression follows from Eq . fl4"7|) by simple substitutions. In- 
deed, the amplitude A A )^ Bl must take the same form of A AB ° Bl when expressed in 
terms of the primed Sudakov variables and primed vectors p'±, defined by 

2 2 
PA' =P- + P + , PB> = P+ + P - , 

s s 
Pi = 0iP_ + a'ip + + p i± > , (51) 
where J_' means transversal to the (pa',Pb>) plane. Since 

2 2 

m A - ql 

p A ' =P- + p+ - q± , 

s 

2 2 

p B ' = p + + b % _ + q ± , (52) 
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we get the following connection between the basis vectors: 



<? 2 

P+ -P+- —P- + q± , 
s 

p_~p_ p + - q ± , 

s 



(53) 



and between the Sudakov variables: 

,?JUPi± 



ft ~ A - ^ + 2 



a- ~ Oj A - 2 

s s 



Pi±> - Pi± + (A - oci)q± + 2 



P+-P- 



(fii -ai)q 2 ± + q±p 



'i± 



(54) 



Keeping only the leading terms, in the kinematical region defined by Eqs. (|3~6"D, 
( |3~7| ) and (|40|), we find that the amplitude A^}^, Bl can be obtained from the RHS 
of Eq.fl4~7j) by the substitution 



k± — > + /3(?j 



<?2± = PBi± -> <?2_L - <?_L 



(55) 



besides the obvious change A —>■ A' and B ^ B' . 



We are now able to calculate the discontinuity ( p2]) . Using Eqs. fl4"T|) and ( |50| ) 
and performing the summation over colour and spin states of the intermediate 
particles, we obtain 



EaAiGBi a *AiC 
•^AB ^A'B' 



GBx 



a u g,bi qL(<i2-q)i 



{A \T C T ] T 3 T c \A)R a 



+fi'i'Jijc{A'\T i T\A)R na + if ij , c (X\T i T*T c \A)R int 



{B'\T j T j \B) 



(56) 



-if ijc (A \T C T> T\A)R int 

Here the first two terms in the square brackets come respectively from the product 
of the abelian and non abelian parts of the amplitudes (|47D and (0); one has 

R a = u(PA>)tf+ [^(k± + (3q 2 ±) - L^(kj_ + pq ± )] + m A 
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x [Lp(k±) - L^(k ± + f3q 2 ±)] j/ + u(p A ) 



(57) 



with 



^(fc±) 



m\f3 2 — k\ 



(58) 



and 



Rna = u(p{ai)$+ [lf(k± + (3q 2 ±) - L»{q l± - (1 - P)q ± )] (j/ Al + m A 



) 



X [-M?1J_) - ^ M (^± + /3?2±)] jt+u(PA) ■ 



(59) 



The remaining terms come from the interference of the abelian and non abelian 
parts. It is not necessary to write them explicitely, because they do not contribute 
to the discontinuity which we are interested in. Remind that we need to calculate 
the helicity conserving part of the discontinuity for the amplitude with colour octet 
state in the t channel and negative signature. As it was discussed just before writing 
the relations (0), we must antisymmetrize the colour factor (S'|T J T^\B) with 
respect to the inversion of the order of the group generators (antisymmetrization of 
the colour factor for the particle A leads to the same result). 

Bearing this in mind, we first calculate the colour factors in Eq . (|56|) . As for the 
R na term, using the relations ([18]) and taking into account that the second of them 
gives us (in the case of adjoint representation {T l )i m = —if 1 j m ) 



ijkjilmjjmn „ Jknl > 




(60) 



we get 



-{B'\T k \B)if kfj f l>f J l3C {A'\TT\A) 

j(B , \T*\B)if^ i (A'\TT i \A) = 
N 2 

- — (B'\T k \B)(A'\T k \A) . 



(61) 
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The same relations allow to show that the colour factors for the interference terms 
are zero. For example, the colour factor of the Ri nt term is 

^if^ j (B , \T k \B)if^ e (A , \T^T°\A) = 

lihfjiB'lT^if,^ {if ljc (A\TT l \A) + (A>\TT c Ti\A)) = 
\{B'\T k \B) ^f klt (A'\TT l \A) + ^if k ,.(A\T^T^ = . (62) 
The same result is obtained for R' int . For the colour factor of the R a term we find 

l -f kfj (B'\T k \B)(A'\ T c T j 'T j T c | A) = 

N i 1 / 

B'\T h \B)(A \T C T T C \A) = -(B'\T k \B)(A \T \A) . (63) 



4 

To obtain this result we used the relation 



T a T c T a = i C2 _ K\t c . (64) 



2 

Here C*2 is the eigenvalue of the Casimir operator T a T a whose expression for quarks 
is 

N 2 -l 



G 



2 



2N 

In turn eq. (|64"D can be obtained using the relations (|i~8|). 

It remains now to calculate R a and R na . It is clear from Eqs.(^) and ( |5P| ) that 
both quantities can be obtained from the general expression 

G(r ± , r' ± ) = u(p A ')^ + I^(r±) (rf Al + m A ) L fl (r ± + r' ± )y( + u(p A ) , (65) 

where r± is linear in k± and r' ± does not depend on k±. Taking into account that 

P+l± = , jf+=0 , Zp+PAi = s(l - (3) , 
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and putting 



we find that the expression 



d(l) 



1 



m 2 A f3 2 - I 2 ' 
may be written as 



(66) 



G(r ± , r' ± ) = s(l - P)d(r±)d(r ± + r' ± )u( PA >) 

~(-m A {3 2 - p/ ± ) 7^ + 2r^] [ 7±/1 (m A /3 2 - /?(/+ f )±) + 2(r + r')±J #>(pa) = 

s(l - f3)d(r ± )d(r ± + r' ± )u(p A ,) [(D - 2) (-m 2 ^ 4 + /3 2 r ± (r + r') ± + m A /?V'±) 

+4r ± (r + r')±(l - /?) - 2m^7'J j/ + u(p A ) ■ (67) 

We have used here the fact that the term f±/'±, being under integration over 
azimutal angles, is equivalent to r±r' ± . 

The term f ± , after the integration is performed, becomes proportional to g[ ± , 
therefore it can be omitted. In fact, because of the relation (see, for instance, 
Eq.(37) in Ref. Q) 



u(p 



u(p A ) = u(p A >){rf A - rf A/ )rf + u(pj 



u(p A ,)(-2mrf + + s)u(p A ) = s (u(p A i)u(p A ) - 2m5x A ,\ A ,) = -isy/^t8x A ,_ Aa , , (68) 

it does not contribute to the helicity conserving part which we are interested in. 
After some algebra, for this part we obtain 



G (+) (ri, r'J = s 2 (l - /3)d{r ± )d{r ± + r' ± ) 

1 1 



\{l-P) + {D-2)^ 



r'f + 



+ Am\P\l-P) 



(69) 



d(r± + r' ± ) d(r±) i 
where the superscript (+) stands for helicity conserving part. 

Let us denote by A^ 3j4 ) the contribution to defined by Eqs. (|TT|) and flT3|), 

which comes from the fragmentation region of the particles A and A'. Applying 



19 



the unitarity condition (^) and Eq.(^) with the colour factors given by formulas 
(0)-(H), we have 



A (3A) = f d ®3(PA+PB]PA 1 ,k,p Bl ) 

s J 27T0? , (q 2 - q) 2 . 



(70) 



where the superscript (+) means, as before, helicity conserving part. The integra- 
tion runs over the fragmentation region of the particles A and A', where the phase 
space element takes the form 

d<S> 3 (p A +p B ;p Al ,k,p Bl ) 1 dp d ( - D - 2 kiLd {D ~ 2) q2±_ 



Here 



2ty 



qi=PA~ PA, 



4s P(l-P) (27r) 2 ^- 1 ) 



(71) 



q 2 = Pb 1 ~ Pb 



Po<P<1 



k = qi - q 2 



Po being a not specified, artificially introduced boundary of the fragmentation region 
of the particles A and A'. It is convenient to split A^ 3 ^ into two parts, 



A (3A) = A (3A) + A (3A) 



(72) 



which respectively contain the terms and R^J ■ Using the phase space element 
(|7l"l), the expressions ([57]) and (l59|) for R^~' and R^/ respectively and Eqs.(|65|), (|69|) 
for G(r±, r' ± ), from Eq . (|70|) (here and below all vectors are (D — 2)-dimensional and 
orthogonal to the (pa,Pb) plane) we get 



A (3A) = ti f d^q, d^q 2 



8 J (27T) - 1 (2tt)^ Jp P q 2 {q 2 - q) 



dp p 2 



/ 2(l-/3) + ^p-2) 



xd(k)d(k + pq)d(k + pq 2 ) 



4 



(g2 ~ qf 

d(k + pq 2 ) ' d(k + pq) ' d(k) 



+ 



+4m 2 A (l - P) {d{k + pq 2 ) - d{k)) {d{k + pq 2 ) - d{k + pq)) 



(73) 
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and 



A (3A) = 9^NH r g^Og dl D -Q q2 f 
na ~ 8 J (2n) D - 1 (2n) D - 1 J 0O 



(2tt) d -i (2tt)^ 7/3 JWq^W 



2(l-/3) + ^(L>-2) 



xd{ qi )d{k + f3q 2 )d{q 1 - (1 - 
4m2,/3 2 



^ 2 



(?2 - g) : 



+ d( 9l - (1 - 

1 ^ + M - d( gi )) (d(k + (3q 2 ) - d{q x - (1 - /%)) 



1-/3 



(74) 



Eqs.(|72|)-(|T4"|) give us the contribution of the fragmentation region of the particles 
A and A'. It should be supplemented by two contributions. One of them comes from 
the fragmentation region of the particles B and B' which, in analogy to Eq. (|72|) , 
can be divided into two parts accordingly to the colour factors, 



A (3B) = A (3B) + A (3B) 



(75) 



Evidently, A^ 3B ) and A^ 3 a B ) can be respectively obtained from Eqs.(|73|) and (|74]) by 
the substitutions 



9l,2 ~> -92,1 , 

/3 — > a = q;g 



q^ -q 
-k 2 



s/3 



/3 — > a . 



(76) 



The other contribution, that we call A^ 3/ \ comes from the multi-Regge region, 
defined by Eqs.([36|) and (|39|) , which is intermediate between the two fragmentation 
regions. This contribution was calculated || and reads 



A ( 



(3J) _ g^NH f d( D -V qi d^ D -^g 2 rP° dp 



x 2 



(2n) D - 1 (2n) D - 1 

q 2 q\ 



(q2~q) 2 



qf(qi — q) 2 q 2 k 2 (<zi ~~ q) 2 k 2 ^ 

This result can be obtained from the expression (|74| ) for A^) by changing the limits 
of integration over (3 and going to small f3 in the integrand. Notice that the abelian 
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part ( |73|) does not contribute in this case. It means that the integrands of Eqs. (|73|) 
and (|74|) are valid in a region wider than that of the fragmentation one, namely not 
simply for f3 ~ 1 but for \t\/s <C /3 < 1. This is not unusual because \t\/s is not an 
artificial bound but a natural one where the applicability should be broken. That 
allows us not to consider the multi-Regge region separately, but to include it in any 
of the two fragmentation regions which, then, become overlapping. Hence, the total 
three particle discontinuity can be obtained as the sum of the contributions of the 
two fragmentation regions, 

A (3) = A (3A) + A (3B) (7g) 

Here the first term in the RHS is given by Eqs.(|72D-(^4D, the other one follows from 
it by the substitution (|76|) , with the parameters (5q and «o satisfying the condition 
s/3o«o = —k 2 . Choosing 

Po = «o = J^f- , (79) 
we get A( 3B ^ from A^ 3A ^ simply by the substitution 

As for the abelian part A^ 3j4 ) of the last discontinuity, it is clear from Eq.(|73[) 
that Po could be put equal to zero because the small (3 region does not contribute. 
This is in accordance with our experience in QED, where it is well known that cones 
of photon emission by two scattered particles do not overlap at high energies. 

Performing the integration over qi in Eq.(|75|) with the help of the Feynman 
parametrization, 

_ Wr+ ^_^^_, (80) 

where 

c(m,l) = [m 2 -/ 2 x(l -x)f-% , (81) 
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we find 



A(3A) 



^r(3-f) 



4 (4 



d/3 



o /3( 5 - D )7o 



7T 2 



92 



M^gl (g 2 - g)< 



^ 2 



2(1-/3) + yP" 2) 



x 



-q 2 + 2g 2 



c(m A ,q) c(m A ,q2) / 



+ 4m 2 A (l-P) 



1 1 

+ 



c(m A ,0) c(m A ,q) c(m A ,q 2 ) 



(82) 



After this step the integration over /3 can be easily performed and yields 

«*r(2-#) 



A(3A) 

^ — */l 



4 (4 



x / dx 
Jo 



7T 2 



D-A 1 
+ 



rf( p - 2 V 
(27r) £ '- 1 g' 2 (g' - g) 2 

2q' 2 



g 



2m^ 



4 D-3J \c(m A ,q) c(m A ,q') 
2 11 



83) 



D-3\c(m A ,q') c(m A ,q) c(m A , 0) / 
In order to calculate the non abelian contribution given by Eq.fl?4D with 



Po — y~k 2 /s, we apply the following trick: we subtract and add this contribution 
considered for a massless quark, A^\m A = 0). While performing the integration 
over qx of the difference 



C> = A£*>-A£?>(m A = 0), 

we may put /3 = 0. With the help of Eq.(^) we arrive at 
fl 4 iV 2 tr(3-f) /• ^~ 2 )g 2 



34) 



5 (3A) 



(47r)f •/ (27r)^g 2 (g 2 - g)^o 
x ^{(2(1-0 + [^(^^ 



/3 



1 



/?)) c(0, g (l -/?)), 



-2g 2 2 



c(m A Ag 2 (l-/?)) c(0, g 2 (l - /?)) 
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4m 2 A {3 2 



1 2 1 

c{m A (3,q{l- (3)) c(m A f3, g 2 (l - /?)) c(m A /5,0) 



The total contribution of the three particle intermediate state to the discontinuity 
is then given by the sum 

A? = Ai 3A) + € A) + A? fl > + € B) + Affi*> (m A = 0) + Ag? (m fl = 0) . (86) 

Here A^ 3 ^ and <^ 3 ^ are given by Eqs. fl83|) and (H) respectively, A< 3B > and tfjgf > 
follow from these equations by the substitution m A m B and A^^m^ = 0), 



which is equal to A^(m B = 0), comes from Eq. (|74|) with (3 = y—k 2 /s. 

5. TWO-LOOP CORRECTION TO THE GLUON TRAJECTORY 

We are now ready to discuss the formula ( |T2| ) for the correction uj( 2 \t). It 
contains the discontinuity As which has been calculated in this paper and, according 
to Eq.([U|), is given by the sum of the contributions A^ and A^ , respectively 
quoted in Eqs. (|2l"l) and (|86|) . Besides that, Eq. (|T2|) involves the leading contribution 
to the trajectory u}^'(t) shown in Eq.((|) and the one-loop corrections to the QQR 
vertex rjjj, which are given by Eqs. ©-(11), (0), (||) and ([31]). Notice that 
the discontinuity A^ itself is expressed in terms of u^(t) and the QQR vertex. 
Therefore, the correction u^(t) looks rather complicated. Fortunately, a series of 
remarkable cancellations occurs among various terms in it. 

First of all, the term A^ 3 " 4 ) in A® cancels all the terms in Eq.(|l^) which contain 
the contribution aq{r, m A ) (see Eq.(p4])) to the QQR vertex T^j defined in Eq . (p2|) . 
It is worthwhile to note here that this vertex enters Eq . fll2|) for uj@'(t) not only 
explicitely, but also through the two-particle discontinuity A^ (see Eq . (|T3|) ) . One 
might easily see the cancellation putting the expressions (p4[) for aQ(r,m A ) and 
(§ for c; (1) (t) into Eqs.© and (0), and comparing the result with Eq.fl53|). Of 

24 



course, the same cancellation exists between the term A^ 3B ) and the terms which 
contain ag(r, m 2 B ). Quite analogously, the term (]85D in Ag cancels all the 
terms in Eq.([l2]) which contain 5 g (r,m 2 A ) fl3lf). In order to make this cancellation 
more evident, one only needs to change the variables in Eq . (|3TD : 

x 2 ^{l-(3)x. (87) 

In turn also the term cancels the terms with 5 g (r, m 2 B ). 

After these cancellations are performed, one realizes that the correction u/ 2 )(t) of 
Eq.(|12"D does not depend on the masses of the scattered quarks. This independence 
was expected because the trajectory, by definition, cannot depend on the properties 
of the scattered particles. In fact, Eq. (|12"D , because of Eqs.(|86D, fl2T|), (|22|) and 
gives us 

g 2 Nt r S D ~ 2 ki 



J 2 \t) = 2At A) (m A = 0) + 



(2tt)^ J q \ ( qi - q f 



10 



(1) te 2 )ln 



(t)Y\n( — )-2^ 1 \t) [a f (t) + a g (t, 



+2a f (q 2 1 )+2a g (qj,0) x ; , f 

where a f (t), a g {t,0) and A^ are given by Eqs.(|), (||), fl2|) and @ 

respectively. 

Looking at Eq. (|88D , we observe that not all the cancellations are performed up 
to this point. Indeed, the trajectory cannot depend on s, therefore the terms with 
ln(s) must cancel each other. Let us note that the cancellation of these terms is a 
consequence of the gluon Reggeization in LLA, which was previously proved 0, 0], 
and can serve as a check of our calculations, whereas the cancellations discussed 
above confirm the gluon Reggeization beyond LLA. 

In order to explicitely demonstrate the cancellation of the terms with ln(s) and 
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to present the correction uj^{t) in a more transparent form, we integrate over (3 in 



Eq. (|74|) with rriA = 0. We go through this step in the following way. We divide 
the integration region into two parts, one from /3 - 



—k 2 /s to 5, the other one 
from 5 to 1, with 8 -C 1 (at the end of the calculation we let 5 go to zero). In 
the first region we may put (3 = everywhere besides the factor d(3/(3; thus the 
corresponding contribution to A^'(m,A = 0) is given by 
g 4 NH fS D -V qi S D - 2 )q 2 1 



8 J (2vr)^- 1 (27r)^ 1 g2 (g2 _ g ^ 



92(9i - <?) 2 (9i 



92 J 



In 



s8 2 



-(9l - 92y 



39) 



In the second region the change of the variable qi 
dependence, yielding 

J ( 1 |(l- /J )«(2(l- /J ) + f(I,-2))=2 ln (i 



^1 — /3)gi factorizes the (3- 



-2ij>(l)-2ij){D-3)- 



2(D-3) 
(90) 



To obtain this result we have used the decomposition 

- 1 dp 



6 (3 



1 f (/(/?) - /(o)) + / X f/(o) 



and the integral 

Given Eqs.(gjJ) 

&gf>(m A = 0) 



D-4 



- 1 



'0 y 

and (|90|), from Eq.fl7|D we get 

g 4 NH r S D ~ 2 ki rf (D - 2) 92 



^(l)-V>(D-3) 



9' 



(91) 



(92) 



8 7 (27^-! (2vr)^ 1 g2( g2 _ g )^ 



In 



92y 



2^(1) - 2tp(D - 3) 



9l(9i - 9) 2 
3 

" 2(£>-3) 



(9i - 92) 



(93) 



We now substitute this expression, together with the expressions (§) for LO^(t), fl23D 



and (|8|) for a/(i) and a g (t,0) respectively, in Eq. fl88|) . Making use of the equality 

D\ -p2 ( D i 

(94) 



S D ~ 2 h 



4rf 2 -f r 2 f-i 



(2 7 r) D - 1 r 2 (r - r') 2 (4vr)f r(D - 3) (-r' 2 ) 



3-4* 
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in order to represent the contribution of a g (t, 0) in integral form, we obtain, as final 
result, 



J 2 \t) 



g 4 N 2 t f S D -Vqi 



4 J (27r) D -y J (27r)^- 1 g 2 2 



In 



(qi-q) 2 (q 2 -q) 2 \(qi ~ <h) 



In 



q\ 



+ 



(qi + q2-q) 2 \(qi-q) 2 J 

x (2^-3) + ^(3- |) -2^(| -2) -V(i) 



(91 - <?) 2 (<?2 - g) 2 (gi + g2-g) 2 , 



(D-3) \A(D- 1) D-A 4 



r 2 



/ dxxil 

r JO 



x 



N (4vr)# (qi-q) 2 



1 



, , • ( 95 ) 

(m 2 — g 2 a;(l — a;)) 2 2 (m 2 — g 2 s(l — x)) 2 2 

This equation gives the two-loop correction to the gluon trajectory in a closed 
form. The summation is performed over quark flavours and t = q 2 . Remind that 
all vectors here are (D — 2)-dimensional and space-like, r 2 = —f 2 . 

6. Summary 

We have calculated the two-loop correction c</ 2 )(t) to the trajectory of the 
Reggeized gluon in QCD. To find this correction we used the scattering process 
of massive quarks at large energies \fs and fixed momentum transfer \/— t. The 
correction is given by Eq.(|T2"D in terms of the helicity conserving part of s-channel 
discontinuity A5 of the amplitude with colour octet state and negative signature 
in the t channel. Furthermore, it depends on the leading contribution u^(t) to 
the trajectory || (see Eq.@) and the one-loop correction to the helicity conserving 
part of the quark- quark-Reggeon vertex T^{t) (see Eqs. (HD-(HD> (13)' (HI) 
and (H)). 
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The discontinuity Ag, obtained in the two-loop approximation, consists of two 
contributions, Af (see Eq.®) and A^ 3) (see Eqs.(||), (H,© and fl74|)), which 



come from two and three particle intermediate states in the unitary condition. 

We have obtained the final result ( p5[ ) for the two-loop contribution tu^(t) 
to the trajectory as a consequence of a series of remarkable cancellations. These 
cancellations exhibit the independence of u^(t) on any specific property of the 
scattered quarks and confirm the gluon Reggeization beyond the leading logarithmic 
approximation. 

Eq.(^) gives us the correction u^ 2 \t) for the space-time dimension D. Of course, 
our main interest relies in the physical case D = 4. Unfortunately, in this case the 
correction shows up divergences which are both ultraviolet and infrared. Indeed, 
the former ones are not difficult to deal with. They are due to simple first order 
poles and can be removed by expressing the bare gauge coupling constant g in terms 
of the renormalized one in the total expression for the trajectory: 



In the MS scheme one has 



11 »r 2 \ 9 



9 = 9* 



2 



1 1 . , . , 1 



-ln(47r) - -?/>(!) 



D-A 2 v ' 2 



+ 



(96) 

where is the renormalized gauge coupling constant at the normalization point /i. 

The infrared divergences, on the contrary, are much more severe. They could 
not be cancelled inside the trajectory u(t) because the gluon is a colour object, 
whereas we may expect their cancellation for the scattering of colourless objects 
only. So, we should check the cancellation of the infrared divergences when the 
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trajectory would be put into the Bethe-Salpeter type equation for the amplitude 
with vacuum quantum numbers in the t channel. We hope to do that in subsequent 
papers. 
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Figure Captions 



Fig. 1: Two particle contribution to the two- loop s-channel discontinuity. One of 
the amplitudes is in the Born approximation, the other one in the one-loop 
approximation. 

Fig. 2: Feynman diagrams for gluon emission. 
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